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ABSTRACT 

We revisit the issue of mixing matrix renormahzation in theories that include Dirac or 
Majorana fermions. We show how a gauge-variant on-shell renormahzed mixing matrix 
can be related to a manifestly gauge-independent one within a generalized MS scheme of 
renormalization. This scheme- dependent relation is a consequence of the fact that in any 
scheme of renormalization, the gauge-dependent part of the mixing-matrix counterterm 
is ultra-violet safe and has a pure dispersive form. Employing the unitarity properties 
of the theory, we can successfully utilize the afore-mentioned scheme-dependent relation 
to preserve basic global or local symmetries of the bare Lagrangian through the entire 
process of renormalization. As an immediate application of our study, we derive the gauge- 
independent renormalization-group equations of mixing matrices in a minimal extension of 
the Standard Model with isosinglet neutrinos. 
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1 Introduction 



One of the most fundamental properties of the well-estabhshed Standard Model (SM) || 
is its renormalizability 0. Renormalizability endows the SM with enhanced predictive 
power that emanates from the fact that ultraviolet (UV) divergences due to high order 
quantum effects can always be successfully eliminated by a redefinition of a finite number 
of independent kinematic parameters of the theory, such as masses and couplings. The 
predictions of the SM have been tested and vindicated with a satisfactory accuracy at 
high-energy colliders, such as the Large Electron Positron (LEP) collider at CERN and the 
Tevatron collider at Fermilab, as well as in low-energy experiments, e.g. in the recent E821 
experiment at BNL where the muon anomalous magnetic moment is measured 0. 

In addition to masses and couplings of the SM particles, however, the quark-mixing 
matrix, the so-called Cabbibo-Kobayashi-Maskawa (CKM) matrix ^ V, needs be renor- 
malized as well [|, |^. In this context, one of the renormalization schemes, most frequently 
adopted in the literature, is the on-shell (OS) scheme of renormalization 0, |], 3, where 
the particle masses are renormalized so as to represent the physical masses at the poles of 
the propagators. It was shown in that the complete UV structure of the counterterms 
(CTs) for the CKM matrix V can be entirely expressed in terms of quark wave-function 
renormalizations. Within this framework, a simple approach to renormalizing V in the OS 



scheme was also presented, consistent with the unitarity properties of the theory ||10|| . 

Even though radiative effects due to the renormalization of an off-diagonal CKM 
matrix were found to be undetectably small in the SM p, |ll|, |T2|, |13[, this needs not 
be the case for its minimal renormalizable extensions. In particular, in |jT^ the above 
formalism of mixing-matrix renormalization was extended to theories that include isosinglet 



neutrinos and so admit the presence of lepton-number- violating Majorana masses [|T^. A 
minimal realization of such a theory is the SM with right-handed neutrinos [|16|, As 
we will further discuss in Section 2.3, in this minimal model the charged and neutral 
current interactions of the W and Z bosons to leptons and neutrinos are described by two 
non- unitary mixing matrices B and C [l^, respectively. Most importantly, the radiative 
effects on the light-heavy neutrino mixing angles contained in the B and C matrices were 
computed to be as large as 15% [jl4|, close to present experimental sensitivities. The SM 
with right-handed neutrinos is an appealing scenario which may explain the smallness of 



the observed neutrino masses and adequately address the solar energy deficit problem [|T? 
through neutrino oscillations |19[. Furthermore, this minimal extension of the SM may give 
rise to a number of observable phenomena, such as lepton- flavour and/or lepton-number 
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violation in /i, r and Z-boson decays |23|, or to possible lepton- number violating 

signals at high-energy colliders, e.g. at the CERN Large Hadron Collider (LHC) p5[] . 



It has been noticed recently [26| that in the OS renormalization prescription presented 
in P], the derived CTs for the CKM matrix naively depend on the choice of the gauge- 
fixing parameter ^ in the class of gauges. This fact is not very desirable, as physical 
matrix elements will be gauge dependent after renormalization. To circumvent this prob- 
lem of ^-dependence of the OS renormalized CKM matrix, several alternative schemes of 
renormalization have been suggested in the very recent literature p6| , 0, ^ . As is 
expected, in all the proposed renormalization schemes, the UV-divergent parts of the CTs 
of the CKM matrix are identical to those derived in the MS scheme ^j. Nevertheless, the 
UV-safe parts of the CTs differ from approach to approach by finite dispersive constants. 
Most interestingly, one may observe that even in the originally suggested OS scheme of 0, 
the gauge- dependent part of the CKM-matrix CTs is UV finite and also has a pure disper- 
sive form, thus indicating the existence of a profound relation between gauge dependence 
and scheme dependence in mixing-matrix renormalization. 

In this paper, we revisit the topic of mixing-matrix renormalization of the CKM ma- 
trix V and of the B and C matrices. In particular, we develop a generalized and manifestly 
gauge-invariant MS approach to mixing-matrix renormalization. The developed general- 
ized MS approach provides a very convenient framework to address the problem of gauge 
and scheme dependences in the existing plethora of differently renormalized mixing ma- 
trices. Moreover, we show how our generalized MS scheme can be successfully employed 
to maintain global or local symmetries of the bare Lagrangian after renormalization. Fi- 
nally, with the help of our generalized MS approach, we can derive the gauge-independent 
renormalizat ion-group (RG) equations for mixing matrices. We explicitly demonstrate the 
theoretical advantages of this method by calculating the one-loop RG runnings of the B 
and C matrices in the SM with isosinglet neutrinos. 

The paper is organized as follows: after briefly reviewing the basic formalism of 
mixing matrix renormalization in the OS scheme in Section 2.1, we present in Section 2.2 
our gauge-invariant generalized MS approach to the renormalization of the CKM matrix 
V, and extend it in Section 2.3 to the renormalization of the corresponding B and C 
mixing matrices in the SM with isosinglet neutrinos. In Section 3, we show how our 
generalized MS approach preserves additional global and local symmetries of the theory, 
which are manifested themselves as sum rules involving neutrino masses and the B and C 
matrices. As an immediate application of our considerations, we derive in Section 4 the 
gauge-independent RGEs of the B and C mixing matrices. Finally, our conclusions are 
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summarized in Section 5. 



2 Mixing matrix renormalization 

In this section, we will first recall the basic analytic formulas for the wave-function and mass 
CTs in the OS renormalization scheme within the context of general fermionic theories, such 
as the SM and its natural extension with isosinglet neutrinos. Then, we will revisit the 
problem of gauge dependence of the OS-renormalized CKM matrix in the SM, and discuss 
its connection to scheme dependence within a generalized gauge-invariant MS scheme of 
renormalization. Finally, our discussion will be extended to the renormalization of the 
mixing matrices B and C that parameterize the neutral- and charged-current interactions 
in the SM with singlet neutrinos. 

2.1 OS renormalization scheme 

In a theory with a number Nj of Dirac fermions, the bare kinetic Lagrangian has the 
following generic form: 

-fnZ]l'\M + 6M)Zl^'h. (2.1) 

In the above, we have employed a matrix notation in the space spanned by the Nf fermionic 
fields, i.e. = (/i, /2, . . . , fNf)- As usual, we adhere the superscript '0' to unrenormalized 
quantities. In ( |2.1| ), the Nf x Nf dimensional matrices Z]/'^ and Z]/^ are the wave- function 
renormalizations for the left- and right-handed fermions, respectively. In addition, M°, 
M and 6M are diagonal Nf x Nf dimensional matrices that contain the bare masses, the 
renormalized masses and their respective counter-terms (CTs). 

The most general form of an unrenormalized fj fi transition amplitude allowed by 
hermiticity P] reads 

E,,(^) = ^PL^f.ip') + ^Pr^Up") + Pl^W) + PR^f:(p") ' (2-2) 
supplemented by the constraints 

^W) = Sg(p') = ^?*(P")- (2-3) 
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In the OS scheme of renormahzation, the wave-function and mass CTs are given by 
^4 = ^2 + m,m^^{rn]) + m,^{m]) + m,T.f:{m])), (2.4) 

= ^2! ^2 (^.^i4K') + '^I^SK') + '^.sJ^K) + "^^sj]*K')), (2.5) 



- ( S,?K) + S^'K) ) - ( S^'(m,2) + ll^im]) ) , (2.6) 

- < ( ^:^) + ^{rnf) ) -mA T^{m^ + ^^{ml) ) , (2.7) 



^rn. = ^m,(S,^(m,2) + S,^(m,2) ) + \ [i:^{n^) + S,^*(m,2)), (2.8) 

where S'(P^) = dTi{p'^) / dp^ and 6Z^'^ are the loop-induced wave-function renormahzations 
defined through the relation: 2']^'^^ = 1 + We should bear in mind that only the 

dispersive parts of the unrenormalized self-energies enter the renormahzation such that 
the hermiticity property of the local Lagrangian is maintained. In the SM, it is S^(p^) = 
mjE^(p^) and S^(p^) = T,j*{p^), and the formulae (p.4|) - (p.8|) reduce to those given in P, P]. 
However, we should stress again that these relations are very specific to the SM and do no 
longer apply to extended theories. 

One well-motivated extension of the SM is the one in which the SM field content 
is augmented by right-handed (isosinglet) neutrinos, thereby admitting the presence of 
Majorana masses in the Lagrangian ^ In this case, the fermionic fields satisfy 
the Majorana constraints: /° = (Z^)*" and = (/i?)*", where the superscript C indicates 
charge conjugation. As a consequence of the Majorana constraints, we obtain the equalities: 

Zl^' = Z]{'\ = ^fnv") ^ S,^(/) = S>^), (2.9) 



where we made the identification S^(p^) = T,fj{p'^). Substituting (|2.9|) into (|2.4| )-(pl 
yields the corresponding wave-function and mass CTs for Majorana fields 



The issue of mixing-matrix renormahzation arises whenever one has to deal with the 
renormahzation of a non-trivial rotation matrix that occurs in interactions relating flavour 
to mass eigenstates. To study this problem, we shall adopt a perturbative framework 



*Here, we have used the symmetry property of the Lagrangian (^]^) under the rephasings, Z^{'^ — )■ 
e'^'Z]^0 and Z^^^ —> e*^'Z]^{^, in order to cast SZ^^ and 5Z^ into a symmetric but fully equivalent form 
than the one presented in [H. 
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in which the classical tree-level Ward identities (Wis) are maintained after quantization. 
As such, one may consider the Background Field Method (BFM) |3y] or the Pinch 
Technique (PT) 



M or even possible diagrammatic generalizations of the latter, 



i.e. the Generalized Pinch Technique (GPT) |^ 



2.2 Renormalization of the CKM matrix in the SM 

As a prototype example, let us consider the charged- current interaction in the quark sector 
of the SM. Specifically, we will revisit the renormalization of the CKM matrix elements Vud 
that enter the vertex transition W~^{p)d{pci) u{pu)- Later on, we will generalize our 
results to the aforementioned SU(2)L(8)U(l)y model with Majorana neutrinos. Within the 



perturbative approaches mentioned above, the following tree-like WI is satisfied [^, [TT|, |26 

p^T^^^'''%p,Pu,Pci) + M^T''^^''%p,p^,p,) 

= -^(^V^,^s:^;\pu)Pl - V^,,PnS,,'/{p,)),{^-W) 

where the summation convention over repeated quark-family indices is implied. In addition, 
in (|2.1CI|) we have defined 



rr"''°(p.p-p^) = + r5:"^'^(P,Pn,P.), (2.11) 

r«"-'^'°(p,p„,p,) = rf'^'-' + Tf'''ip,Pu,P,), (2.12) 

Sd'YiPd) = - m', + , (2.13) 

where Pjj^ and are the bare W~^ud- and G^ud- couplings at the tree level, 

and F]^^"'^(g,p„,prf) and Ff^"'^(g,p„,prf) are the corresponding higher-order unrenormal- 
ized one-particle irreducible vertices evaluated within e.g. the PT or the BFM. Similar 
identifications also apply for the unrenormalized two-point correlation functions S^^i{pu) 
and S%d{pd). 



Following the procedure outlined in ||3^, we require that the same tree-level WI ( p.lO|) 



which involves unrenormalized quantities only, holds exactly true after renormalization. 
This condition can be successfully enforced within the gauge-independent MS scheme of 
renormalization. Nevertheless, in any other favourable scheme of renormalization, the 
renormalized parameters of the theory will differ from those in the MS scheme by UV-finite 
constants. The renormalized quantities may be determined in terms of the unrenormalized 
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ones through the relations: 
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= + 5 Ml, , 

^0 = V + 6V, (2.14) 

[P,Pu,Pd) — ^Luu' ^Ld'd^ fj. [P,Pu,Pd), 



^^^\p,p^,p,) = Z'j:zlO^Z'J!r^^^''^''%p,p^,p,), (2.15) 
S^lM) = Zl!SS^:'{p,)Zl!' (withg = n,rf), (2.16) 



where Z^l^. = Z\I^.^^Pl + Z\l^_^^Pii and 5V stands for the mixing- matrix renormahzation of 
the CKM matrix V . In ( p. 15 ), we required that W^ud- and G^ud- couphngs be UV finite 
after the external wave-function CTs for the ly-boson, the would-be Goldstone boson G"*", 
and the u- and ci-type quarks have been properly taken into account. 

Since our main interest is to compute the UV divergent part of the W^ud vertex 
in the presence of flavour mixing and so determine the UV-divergent structure of the 
mixing-matrix CT 5V ^ we shall therefore focus our attention only on the chirally-projected 
WI ( |2.10| ) related to the expression PRFj^^^'^'^Pi. In particular, we have 



^1/2 ^«,l/2t ^ ^>.^W+nd,0 ^ ^G+Ud,0 ^ p^^d,l/2 



= -Zll'z,^^zf/'^ Pn{s~''\p^)V' - V S-'\p,)) P,zt"\ (2.17) 

To simplify notation in ( p^.l7D , we have employed the matrix representation for the quark 
wave-functions and their inverse propagators, i.e. Zf^^^"^ = Z]l^^^^ and S^^'^{pq) = S~^'^^{pq), 
with q = u,d. Substituting (|2.15|) and ( p. 161) into (|2.17|) gives 

J- M + ^^^W J- ) - - 

xPn(S-\p^)Zl'-'/'v'zi'^' - Zf /'W^zt'/'^ S-\p,)) Pl. (2.18) 



The requirement now that the WI (|2.18|) retains its original form (|2.10| ) where all quantities 
are replaced by their renormalized ones gives rise to the following consistency conditions: 

A/2 _ ry-\ ry _ 7 , 



Z'^' = Za+ = Zw[l + ^), (2.19) 

V = Z2'~'/V°Z^''/' = ^^'V2ty0^d,-i/2t_ (2.20) 

The two equalities in ( p.l9| ) are exactly satisfied within the PT and BFM frameworks 

26| . The double equality in (|2.2CI|) assures the unitarity property of the renormalized mixing 
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matrix V, i.e. V ^ = V^. Most importantly, ( |2.2(j| ) determines the analytic structure of 
the CT 6V. Employing the usual decomposition for the wave-function renormalizations, 



i.e. Z^' =1 + 1 and Z^' ' = 1 + | 6Z'^'^, we arrive at the perturbative and more 

familiar form for 6V 

6V = ^{SZ""'^ - - ^V(^6Z'^'^ - 6Z'^'^^), (2.21) 



where the u- and d- quark wave-functions have to satisfy the constraining relation ITT], ^ 



+ 6Z'''^^)V = ^V[6Z'^'^ + 5Z'^^^^). (2.22) 

In the absence of flavour mixing, i.e. for V" = 1, this last relation simplifies to the known 
one Z"'^ = Z"'-^ [0, g . 

Several important remarks and observations regarding mixing-matrix renormalization 
are now in order: 

(i) The UV poles of 5V are entirely specified by the wave-function CTs of the left-handed 
u- and d- quarks to all orders in perturbation theory. Moreover, with the definition 
of 5V in ( p.21|) , V is automatically unitary through the order considered. 



ii) The LHS of the WI ( |2.18D is gauge-independent, when the chirally-projected ampli- 



tudes PFjr^^'^'^{p,pu,Pd)PL and PRr'^^'^'^{p,pu,Pd)PL are evaluated by setting the 
external particles on their mass shells. Consequently, the RHS of ( p.l8| ) must 
be gauge-independent as well. This can only happen, if V and hence 6V are 
gauge-independent |jlT], |26|. For example, unlike in the MS scheme 0, 6V is ^- 



dependent [p!2| , \L3i , |27| , Pq] in the OS scheme of renormalization. As we will see below 
in (iii), however, because the gauge- dependent part of SV is UV finite and has a pure 
dispersive form, the ^-dependent terms of an OS renormalized CKM matrix V can al- 
ways be related to finite gauge-independent constants in a generalized and manifestly 
gauge-invariant MS scheme of renormalization. 

(iii) There exists an underlying symmetry in the renormalization of V^, reflecting the 
presence of a general intrinsic freedom in redefining mixing matrices at higher orders.^] 



^Here, we should remark that the analytic results for the wave- function CTs and 6V obtained within 
the (G)PT or BFM are identical to those derived within the conventional framework of i?^ gauges [l^ , 
with the additional restriction that the gauge-fixing parameters related to the photon and the Z boson are 
equal, i.e. = ^.z- 

^The existence of such a scheme dependence in the renormalization of was first pointed out in fl^ . 
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The presence of this higher-order scheme arbitrariness in the renormahzation of 
may be described as follows. We know that the CKM matrix is the product of 
two unitary matrices [/]*'° and f/^'° relating the weak to mass eigenstates of the left- 
handed u- and d-tjpe quark fields, respectively, i.e. V'^ = [/]^'°[/^'°^ If we now perform 
the following perturbative shifts in the left-handed quark wave functions and their 
mixing-matrix CTs: 

SZ""'^ SZ""'^ + c", SZ'^'^ SZ'^'^ + c^ (2.23) 

where and c"^ are anti-hermitian, gauge-independent UV-finite constant matri- 
ces, i.e. c"''^ = — c"''^'!" and c"''^ = 0{SZ'^'^,5Z'^'^), then none of the important key 
equalities (p.20|) , ( p.21| ) and ( |2.22| ) will change through the order considered. On 



the basis of the above unitarity symmetry, we may generally define the manifestly 
gauge-independent CKM matrix CT: 

SV = 6V^ - ^c"\/ + ^\/c^ (2.25) 

where 6V^^ is the corresponding CT in the gauge-invariant MS scheme evaluated 
in ( |2.21| ). Evidently, 6V does not give rise to gauge dependences in the computa- 



tion of physical transition amplitudes, e.g. in the top-decay amplitude t ^ W + b. 
However, the UV-finite constants and c'^ introduce a scheme dependence into the 
renormahzation of the CKM matrix Thus, all the various schemes of renor- 



mahzation proposed in the literature |2g, 0, |13|, including the OS 

scheme may be represented by the gauge-invariant expression ( |2.25| ) with 

appropriate choices for the UV-finite matrices c" and c'^. 

Finally, it is interesting to notice that the number of independent parameters con- 
tained in the anti-hermitian matrices and c'^ is exactly equal to the number of 
group parameters, the so-called mixing angles, that generate the unitary rotations in 
the left-handed u- and d- quark flavour space. To elucidate the above point, let us 
consider the renormalized unitary matrix in the left-handed w-quark flavour space, 
U^. In particular, can be written as 

Uliei) = exp (z^^T"), (2.26) 

where 6'^ are the group parameters and are the associate generators of the U{n) 
flavour group, satisfying the usual Lie-algebra relations: [T^,T^] = if"-''^T^, with 
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0. If we now shift 9'^ by a finite higher-order amount 



59'^, then the unitary matrix ( p.26| ) exhibits the following variation: 

1 



be 



6 



jbdx jcexQdQe 



" ebcdnd 
2/ 



69tr 



Ul{9l) + 0{59l') . (2.27) 



In writing the RHS of ( |2.27| ), we have employed the Baker-Hausdorff formula for 
infinitesimal non-Abelian rotations. Comparing (p.27|) with ( p. 241) , we immediately 
recognize that the anti-hermitian matrix is given by 
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jacx jbdxQC Qd _j_ 



arjib 



591T 



(2.28) 



where the ^^°'s parameterize the scheme-dependent shifts in the mixing angles 9'^. 
Notice that the anti-hermitian matrix c" in (|2.28| ), determined by means of (|2.27|) , 
preserves the unitarity properties of by construction. Analogous determinations 
may also be found for and hence for the scheme-dependent part of V in (|2.25). 



2.3 Mixing renormalization in the SM with singlet neutrinos 



We will now discuss the problem of mixing-matrix renormalization in an SU(2)(8>U(l)y the- 
ory with a number Nq of fermionic doublets and a number Nr of right-handed (isosinglet) 
neutrinos. The interaction Lagrangians of this model which describe the couplings of the 
W^, Z and Higgs {H) bosons to the Nq charged leptons, /j, and {Nq + Njij Majorana 
neutrinos, n^, are given by 



9u 



4c^ 

9^ 



Z,n,Y(a,PL - C*^PR)n, 



(2.29) 
(2.30) 



4M, 



Hui 



w 



[midj + mjC*^)PL + {miC*j + m^Cj- ) P^] n^- . (2.31) 



Here, we follow the conventions of JT^, |2g. In (p:29|)-(|23l|), B and C are Nq x {Ng + Nr) 
and {Ng + Nji) x [Nq + A'^R)-dimensional mixing matrices, defined as 



B, 



Ng 

k=\ 

Ng 

k=l 



(2.32) 
(2.33) 
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In (|2.32D, the Nq x Nq unitary matrix V occurs in the diagonahzation of the charged-lepton 



mass matrix and relates the weak to the mass eigenstates of the left-handed charged leptons. 
Correspondingly, the [Nq + Nji) x [Nq + Nn) unitary matrix f/" in ( p.33|) diagonalizes the 
symmetric neutrino-mass matrix 

M" = f "° ) . (2.34) 

y niM J 

through the unitary transformation 

jjnT^njjn = j^n = diag(mi, mg, niNa+Nn) ■ (2-35) 

At this point, we should recall again PI that the mixing matrices B and C obey a number 



of basic identities which ensure the renormalizability of the theory: 



Ng+Nr 

Bi^B*, = 5u', (2.36) 

k=l 

Ng+Nr 

CikCkj = ^ij ^ (2.37) 

k=l 
Ng+Nr 

X! BikCki = Bli, (2.38) 

k=l 

Ng 

Y.B*uBi, = a,, (2.39) 
1=1 

Ng+Nr 

mkBikBvk = 0, (2.40) 

k=l 
Ng+Nr 

Y mkBikQk = 0, (2.41) 

k=l 

Ng+Nr 

E ^kCikC.k = 0. (2.42) 

k=l 



The last three relations (|2.40|) -( p.42|) are manifestations of the presence of lepton-number 



violation in the neutrino sector. Instead, if theory conserves lepton number, these three 
identities are not necessary. In this case, Majorana neutrinos are either degenerate in pairs 
forming massive Dirac neutrinos or unpaired giving rise to massless Majorana- Weyl two- 
component spinors. As a consequence of lepton-number conservation, the mixing matrix 
elements C*- are absent from the ZniUj- and HniUj- couplings in ( |2.3CI| ) and ( p.31|) , so the 
Znjnj-coupling becomes purely chiral, proportional to the operator nij^PLUj. 
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The renormalization of the mixing matrices B and C can be carried out in a way very 
similar to the SM case. Following analogous steps for the one-particle irreducible vertex 
functions T^'^^j and r^"»"^ in the BFM or PT, we find 



B 
C 



yl-\l2 j^O ^",1/2 

L L 
771,-1/2 ^0 2^n,\l2 



U/2t ^0 -V2t 
,ri,l/2t ^0 ^n,-l/2t 



(2.43) 
(2.44) 



where Z^^'^^ and Z^'^''^ are wave-function renormalization matrices for the left-handed 
charged leptons and Majorana neutrinos, respectively. Equations (|2.43|) and ( p.44| ) lead 
perturbatively to the mixing-matrix CTs 



bB 



5Z^'^^ ] B 



1 
4 



B ( 6Z'''^ - 6Z"'^^ 



6Z 



n,L] 



c 



4 



C i 6Z'"'^ - 5Z^'^^ 



In addition, the following constraining relations are satisfied: 



6Z 



+ ) B 



2 
1 



B { + 



C { 5Z''^^ + 5Z 



(2.45) 
(2.46) 

(2.47) 
(2.48) 



It is important to observe that the renormalized mixing matrices B and C given in ( |2.43|) 
and (|2.44|) as well as the perturbative definitions of the mixing-matrix CTs 6B and 5C by 
means of ( |2.45| ) and ( |2.46| ) fully satisfy the identities ( |2.36D -( p.39D . However, the compati- 
bility of 5B and 5C with the remaining identities ( |2.40| )-( p.42| ) proves more subtle and will 
be discussed in Section |^. 

In the MS scheme, the mixing- matrix CTs 5B and 5C become gauge independent, 
only after the corresponding gauge-dependent part of the tadpole graphs are included [^, 
|30| , |27[| . To further illuminate our procedure, we calculate the Higgs-boson tadpole 
graph r-^ induced by the boson, and the associate would-be Goldstone boson and 
ghost fields c"*", c^, i.e. 



rf^)(o) 



9n 



Ml 



327r2 Mw 
. 9w 



167r2 



- l)Mw 



Ml 



B^{S),iwMl,iwMl) 
(l + 5o(0, M^,M^ 



(2.49) 



where the one-loop function BoijP ,m\^m^ is defined in Appendix A. From (|2.49|) it is 
easy to see that only the M^-dependent part of the tadpole depends on and should be 
included in the scalar part of the self-energy transitions (cf. ( p.2|) ). More precisely, the 



12 



M^- dependent part of the tadpole graph effectively induces a gauge-dependent shift to the 
i/-boson VEV v. 

''T = sH;^ = + MO,i«^MU.vMl)). (2.50) 

Similarly, the Z-boson loop causes an analogous gauge- dependent shift to v, i.e. 

V)" - 211^^ - 3£|^-(' + B.(0,(MMI)). (2.51) 

Then, the ^-dependent VEV shifts ( |2.50| ) and ( ^.51| ) contribute the following term to the 
scalar part of the Major ana- neutrino self-energy transitions: 

Ej^'*^^ Pl = - (^)^"'' {m,Q, + m,C:^ ) . (2.52) 

Note that if neutrinos are Dirac particles, one has just to drop the term proportional to 
C*j on the RHS of ( p.52| ). As we will see more explicitly in the next sections, the tadpole 
contribution ( p.52|) plays an instrumental role to render the mixing-matrix CTs 6B and 6C 
gauge independent. 



3 Neutrino mass-mixing sum rules 



As was already mentioned in Section ^ the neutrino mass-mixing sum rules ( |2.4(j| )-( p.42| ) 



are very essential to ensure the renormalizability of the theory. These sum rules are obtained 



by projecting out the zero texture in the Majorana-neutrino mass matrix (|2.34|) as follows: 



Ng+Nr 

Y: m.UllUp, = {U-M-U-^\, = M-r = 0, (3.1) 
fe=i 

for /, /' = 1, 2, . . . , Ng- The zero texture is protected by the gauge symmetry of the theory, 
since the contributing 5-dimensional gauge-invariant operator L^^^^L'-^ is absent from 
the local renormalizable Lagrangian, where L and $ are the lepton and Higgs doublets, 
respectively. This operator is radiatively generated at the one- |]3^, [1^ and two- loop 
levels and is UV finite. The neutrino mass-mixing sum rule ( p.l|) is no longer valid, if 
the theory is extended by one Higgs triplet A/,, since the afore-mentioned lepton- number- 
violating operator can now appear in the tree level Lagrangian through the term L'^A^L'-'. 

In the following, we will show that renormalization of f/" does not spoil the basic 
identity ( p.l| ) in the MS scheme. Within the scheme of renormalization outlined in Section 0, 
the CT matrix 5f/" of f/" is given by 

= if/"(^5Z"'^3^ - (3.2) 



13 



Observe that ( p. 21 ) may also be derived by setting = 1 in ( |2.45| ). In order that the bare 
and renormahzed mixing matrices [/"'° and [/" obey ( |3.1| ), one has to show that 

(5f/"M"f/"^ + U''W'6U''^ + U''6M'^U''^)ii, = 0, (3.3) 

namely the corresponding mixing and mass CTs obey also (|3.1|) up to higher orders of 
perturbation proportional to ((5f/")^. 

Before offering a proof of ( p.3|) for the most general case, let us first gain some insight 
from considering an one-generation model with one right-handed neutrino only. The mass 
spectrum of this simple model, which essentially resembles the known seesaw scenario 



consists of two Majorana neutrinos: a light neutrino u observed in experiment and a yet- 
undetected superheavy one A^. Most interestingly, in this model the elements of the bare 
mixing matrix f/"'° are entirely determined by the two bare mass eigenvalues, m° and m%, 
of M'^'° in dpI): 



rrO _ .Ja^TjO ttO _ TjO ttO _ jjO /o 



where we have dropped the superscript 'n' from f/" and have chosen the phase convention in 
which the elements U^^ and W^j^ are positive. In the MS scheme, the mass, wave-function 
and mixing CTs are found to be 



'^'"^ " levr + m;v V ^ ^""^ 



= - 77- _ , _ 77^- 772 ^TT- ^'uv , (3.5) 





««- ^ / 
lOTT \ 
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O^w ^ I 
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Stt V 


szk. = 


(^w ^ ( 
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lOTT 



167r + m^v V 

\. , . mf ml m^mN \ 

iz ml 3ml m^ \ 

^ 24 + M2,m. + m;vy ''''' 

^ 2c2 + M^^m. + m^vJ ''''' 



J- . 3a^ (mjv - m^) 



2 



In the above, Cjj = UyiU*^ according to ( p. 331 ) (with i,i = u, N), mi is the charged-lepton 
mass, and Cuv is an UV constant defined in Appendix A. As was discussed in Section 0, 
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we find that the mass CTs 6m,^ and 6m n, and the CT matrix SU"' computed by ( ^.21 ) are 
gauge-independent only after the tadpole contributions are included. Moreover, these CTs 
satisfy the basic identity (|3.3|), i.e. 

6m„U^^ + dmNW^N + '^.m^.U^y 5Uyy + 2mNUyN SU^n = 0. (3.8) 



In the above simple Majorana-neutrino model, it is still possible to follow an alterna- 
tive approach. Specifically, the same results would have been obtained if we had considered 
the elements of U"''^ as functions of m° and m^, i.e. f/"'° = f/"'°(m[l, m^). In this case, 
the mixing-matrix CT 5f/" is calculated as 

5f/" = 5rn„ h OTJiN , (3.9) 

only orriN 



and the basic relation (3^) of the CTs will be satisfied by construction, even within the OS 
scheme of renormalization. In addition, 6U"' defined in terms of Sm^ and drriN is gauge- 
independent ||2^. Instead, if we had employed (|3.2| ) to compute 6U"' in the OS scheme. 



the resulting expression would have naively been gauge- dependent and have violated the 
CT relation ( |3.8| ) by UV finite terms.| Nevertheless, even if the mass CTs are evaluated 
in the OS scheme, we can always restore the validity of the sum rule ( ^.1[ ), along with the 
constraining relation ( p.3| ), by redefining 6U"' in a gauge- invariant manner. To be specific, 
exactly as we did in Section |], we add a gauge-independent and UV-finite anti-hermitian 
matrix c" to the MS CTs of f/": 

SU"" = + ^c"f/", (3.10) 

with c" = — c"^ In agreement with our phase conventions for the matrix elements of t/" 
in (|3.4|) , it is sufficient to assume that the matrix takes on the form 




(3.11) 

where c is a real constant. Then, the parameter c can be uniquely determined by requiring 
that the constraining relation ( |3.8| ), with the mass CTs Sm^ and ^m^r computed in the 
OS scheme, holds exactly true. In this way, we were able to verify that the so-derived 
mixing- matrix CTs 6Uyu and 6Uun are identical with those obtained by virtue of (p.9|) . 

For the more realistic case, for which the SM contains more than one right-handed 
neutrino, the unitary matrix [/" cannot be entirely expressed in terms of neutrino masses. 



3 This last result was earlier observed in |l^ 
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and the alternative approach based on ( ^.91 ) turns out to be not very practical. Instead, one 
may utilize the more general approach described above, in which anti-hermitian constants 
are added to the MS CT SU"^. Within this generalized framework of the MS scheme, the 
constraining relations (|3.3| ) reduce the number of independent constants c". The remaining 
freedom should be fixed by comparing the theoretical predictions for observables involving 
the undetermined matrix elements of with experiment. Here, we should stress again 
that the anti-hermitian constants c" are only required if the neutrino-mass renormalizations 
6mi are computed by ( |2.8| ) in the OS scheme. 

Unlike in the OS scheme, in the MS scheme one has to pay the price that the MS- 
renormalized neutrino masses are not the physical ones, namely the poles of the neutrino 
propagators. However, as we will now show, all the basic symmetries of the theory, including 
the one reflected in the sum rule ( |3.1|) , are preserved and the addition of anti-hermitian 
constants c" is no longer needed. In particular, we will provide a general proof of the 
validity of the constraining relation ( p.3| ) governing the neutrino- mass and mixing CTs in 
the MS scheme. To this end, we first substitute f|3.2|) into fB.3|) 



nT ^ ^ 




r 1 


jjn 






.4' 






jjn 


[2 







n i\/rn xrmT 



?n T rnT 



4 



u 



nT 



n ^M*,\JV jjnT 



Jjn 



UV ^ JjnT 



U' 



nT 



(3.12) 



where the superscript UV on and indicates their UV divergent parts, and /, /' 
1,2,..., Ng- In deriving the last step of (p.l2|) , we have used the relations 



^A/,UV 



2 2 



(3.13) 
(3.14) 



These relations may be straightforwardly obtained with the aid of (|2.4| ) and ( p.8|) (see also 
(4.8) and (4.9) in 0). The first term (?7" S*^*'UV ^7^^^^^, ^^le RHS of the last equality 
in ( p.l2|) vanishes by itself, as a consequence of the absence of the operators ^'zl^'Fl from 
the local Lagrangian. The second term vanishes, only if 



^7" M'^S 



0, 



(3.15) 
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or equivalently if 



al' 



0, (3.16) 



with a = No + 1, No + 2, . . . , No + Nr. Equation ( |3.16| ) is derived by inserting the unity, 
jjnTijn* ^ between M" and S^-^v in ( glsD , and noticing that (?7" t/"'^);. = 

Ml^* = m*jji^. Employing the analytic expressions for the neutrino self-energies in Appendix 
A (see also ( [4 .91 ) below), it is not difficult to show that ( p.l6|) is indeed valid. In fact, the 
vanishing of (f/"* S"^'^^ results from the absence of the lepton-number-violating 

kinetic terms i^iLi^^aR dimension 4 from the original Lagrangian in the flavour space. In 
the SM with right-handed neutrinos, the violation of lepton number occurs softly through 
the Majorana operators i^aR^pR of dimension 3, which is reflected in the neutrino mass- 
mixing sum rule ( |3.1| ). This completes our proof of ( |3.3| ). 

We end our discussion by remarking that our general approach to the mixing-matrix 
renormalization may be applied to supersymmetric theories as well, e.g. to the unitary 
matrix that diagonalizes the neutralino mass matrix in the Minimal Supersymmetric 
Standard Model (MSSM). In this case, a convenient framework for renormalization that 



respects supersymmetry is the so-called modified dimensional reduction (DR) scheme |39 



Alternatively, one may work in the MS scheme and translate the results into the DR scheme. 
As in the Majorana-neutrino case, the particular zero texture in the neutralino mass matrix 
is protected in the DR scheme, but needs be reinforced by adding appropriate anti-hermitian 
constants to the DR-renormalized neutralino mixing matrix, if the neutralino masses are 
renormalized in the OS scheme. 



4 Renormalization-group equations 



As an immediate application of our study in Sections | and |^, we derive the gauge- 
independent RGEs of the mixing matrices B and C in the SM with right-handed neu- 



trinos. With this aim, we first compute their respective CTs 6B and 6C by means of ( |2.45|) 
and (l2l46| ) in the MS scheme [|0 



5Br 



Ng+Nr r 

- E 



No R 

^, 2{mf - ml) 



mf + mj,) Sjlf + 2mimii Sjjf + 2 {miT!iiY + miiT}{,\ 



LR 



Ik 



k+i 2(m2-m2) 



(m^ + ) T.^^ + 2mkmi S^^ + 2 ( mfcE^! + miS^^ 



uv 



uv 



(4.1) 
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kj 



m'f + ^k) Kk^ + 2mimfc S^fc"" + 2 ( niiE^'' + m^E 
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^ki 
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Ng+Nr 

- E 



ik 



2(mi — m3) 



{ml + m^) E^f + 2mkmj E^f + 2 ( mfcS^j" + m^Sj^ 



n,D* 



UV 



• (4.2) 



Note that the expressions ([4.1|) and ( [4.2|) pertain to Dirac neutrinos. In case of Majorana 
neutrinos, these expressions are supplemented by the constraints stated in (p.9|), with T,f- 



replaced by S^^ . 

From the analytic results presented in Appendix A, it is straightforward to deduce 
the UV-divergent parts of the self-energy functions occurring in ( [4.1D and ( [4.2| ) in the 
gauge. We start listing the UV-divergent parts of the individual self-energy functions for 
the charged leptons 



K^W ) 



a,, 



IGtt 
+ BuBi, 
a 



5w Asli^ + 2iw + 



2 



UV 



IQrc 



2mf 



UV 



(4.3) 
(4.4) 



167r 



5,, ( 4.2 (3 + e,) - ^^^ii^^ (3 + Cz) +iw + ^ 



+ 2BiiBl 



^1 

* /\/f2 



Cl 



UV 



(4.5) 



Here and in the following, we consider the summation convention over repeated indices in 
their whole allowed range, unless explicitly stated otherwise. Specifically, charged lepton 
indices, such as / and are summed from 1 to No, and neutrino indices, e.g. k, n, from 
1 to Ng + Nfi. The divergent pieces of the self-energy functions for Dirac neutrinos are 
given by 



■.n,L\UV 
-ij ) 

,i?\UV 
ij ) 

n,DNUV 
ij ) 



a,. 
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Cij [ 2^w + % ) + BiiBij 
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yj UliUlj 



uv 
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rrii 
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2cl 



+ '^Bi^Bij 
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uv • 



(4.6) 
(4.7) 

(4.8) 



n,LxUV 
ij ) 



a,, 



If the neutrinos are Majorana particles, the UV parts of the self-energy functions then read 

Cuv, (4.9) 



+ CikCkj 
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Ml 
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/^n.AfxUV 
[^ij ) 



167r 



{midj + mjC*j) + 



+ 2{miBlBij + mjBuB^j 



mf 



w 



Cuv • 



(4.10) 



The above analytic results of the UV pole structure of the neutrino self-energies 
reveal that the RG running of the mixing matrices B and C depends on the nature of 
neutrinos, namely on whether they are Dirac or Majorana particles. In the MS scheme, 
the /i-dependence of B and C may be computed by the beta functions [3b and (3c as 



dB 

Pb = f^-r- 

djji 



lim eg^^ 



^ 5B^\ 



dC 

PC = f^-r 



lim eg^^ 
e^o og. 



(4.11) 



where we have employed the fact that iidg^/dn = —egw + 0{gl^). With the help of (^4.11|) 
and of ([4.1| ) and (|4.2| ), we obtain the following one- loop beta functions for the Dirac neu- 
trinos: 
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(4.12) 
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(4.13) 



Correspondingly, the one-loop beta functions for the Majorana neutrinos are given by 
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a, 



■W 
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(4.15) 



It is worth commenting again on the fact that the beta functions Pb and (3c iii ( |4.12| )- 
( [4.15| ) become gauge independent in the MS scheme, only after the gauge-dependent tadpole 
terms proportional to Mfj have been added to the self-energy functions (or S^^). To 
the best of our knowledge, the beta functions Pb and Pc represent the most general results 
pertaining to the one-loop RG-running of the mixing matrices B and G* in the existing 
literature of the SM with isosinglet neutrinos. However, we should remark that the derived 
RGEs for Pb and f3c are only valid for energies larger than the heaviest neutrino mass. We 



have not considered the threshold effects due to the decoupling |^ of the heavy neutrinos, 
as these effects highly depend on the particular low-energy structure of the model ESI, HI, 
44|, ESl O] and will therefore be studied elsewhere. 



5 Conclusions 



We have revisited the problem of gauge dependence that occurs in the renormalization of 
mixing matrices, within the context of two generic frameworks: (i) the quark sector of the 
SM and (ii) the leptonic sector of the SM with isosinglet neutrinos. Although we confirmed 



the earlier observations p6| , [T^ , p!3| that an on-shell renormalized mixing matrix contains 
gauge-dependent terms, we have observed, however, that these terms are UV finite and have 
a pure dispersive form. Because of this last fact, we have found that these naive gauge- 
dependent terms can always be absorbed into the definition of a manifestly gauge-invariant, 
but physically equivalent, mixing matrix, where the latter is evaluated within a generalized 
MS scheme of renormalization. This generalized scheme of renormalization is obtained by 
adding gauge- independent anti-hermitian constants to a gauge-invariant, MS-renormalized 
mixing matrix (cf. ( |2.25[ ) and ( p.lOD ). In this way, the different renormalization schemes 
proposed in the literature 



II, PI, 13, III, including the OS scheme 
be described for appropriate choices of the anti-hermitian constants. 
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Our generalized MS approach to the renormahzation of mixing matrices may also be 
conveniently applied to maintain fundamental global or local symmetries of the unrenor- 
malized Lagrangian. For instance, our approach may be utilized to protect the texture- zero 
structure of the Majorana- neutrino mass matrix ( 2.34 ) or similar constrained structures of 
predictive neutrino-mass models. Even though such additional symmetries are automati- 
cally preserved in the MS scheme, they become distorted by UV-finite terms and so need be 
reinforced, if the renormalized masses are evaluated within other renormahzation schemes, 
such as the frequently adopted OS scheme. Most importantly, our approach of mixing- 
matrix renormahzation may be applied to supersymmetric theories as well. In this case. 



the corresponding generalized DR approach may be used to renormalize the mixing ma- 



trices that occur in the chargino and neutralino sectors p8|, as well as in the squark and 



Higgs scalar sectors p], |3 of the MSSM. 

As a byproduct of our study, we have derived in Section ^ the gauge-independent 
RGEs for the MS-renormalized mixing matrices in the SM with isosinglet neutrinos. The 
so-derived RGEs are valid for energies that are higher than the mass of the heaviest of 
the heavy neutrinos. We have not taken into account the decoupling effects due to heavy 
neutrino thresholds, since they crucially depend on the low-energy structure of the model. 
However, they prove important to properly describe the RG-running of the observed neu- 
trino masses and the neutrino-oscillation angles at lower energies. We plan to return to 
this phenomenologically interesting topic in the near future. 
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A Neutral and charged lepton self-energies 



Here, we present analytic expressions for the neutrino and charged-lepton self-energies in 
the renormalizable gauge. The Feynman diagrams that contribute to the neutrino self- 
energies are shown in Fig. 0(a)-(d), while the corresponding graphs giving rise to charged- 
lepton self-energies are displayed in Fig. |l|(e)-(g). Our analytic results are expressed in 
terms of the usual Pasarino-Veltman one-loop functions 



where the Minkowski space is extended to D = 4 — 2e dimensions and fi is the so-called 
't-Hooft mass scale. Also, we adopt the frequently-used 4-dimensional convention for the 
Minkowskian metric g^'^: g^'^ = diag(l, — 1, — 1, — 1). The one-loop functions Bq and S^, 
defined in ( |A.1| ), are given by 

Bo{p\ ml, ml) = Cuv - In (^) + 2 + 1 [ (m^ - m?) In 



\ fl"^ J p 
+ X^^'^{p'^,m'l,ml) cosh 



I f ml + m\ — p^ 



2mim2 J J 

B^ip'^,ml,ml) = p^Bi{p^,m\,ml) , (A.3) 
with Cuv = ^ — + ln47r, A(x, y,z) = {x — y — z)"^ — Ayz and 

B^{p\mlml) = {b^{p\ m\, ml) - Eo(0, rnl m^)) - ^ B^{p\ m\, ml) . (A.4) 

The one-loop function B^ij^^ml, ml) evaluated at = simplifies to 

B„(0M,'nl) = C„v - 1„ (!^) + 1 + '4±^ In f^) . (A.5) 

From this last expression, a useful identity relating the arguments of the i?o-function at 
= may easily be derived 



„ n o^ m. 



So(0, mf,m^) = , ' „ Eo(0, mf,mf) So(0, m^, m^) + 1. (A.6) 

— m2 mf — m2 

Equation (|A.6|) may be successfully employed to check the gauge independence of physical 
quantities. 

We first derive analytic expressions for the case of Dirac singlet neutrinos. In the 
gauge, these are given by 

'25i(/,mf,M^) + Boip^mlM^) + 1 
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Figure 1: Feynman graphs contributing to (a)-(d) neutral and (e)-(f) charged lepton self- 
energies in the unitary gauge. If neutrinos are Dirac particles, the graph (b) is absent. 
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Note that in ( [A.9| ) the p^-independent terms represent the ^-dependent part of the tadpole 
contributions. These contributions are crucial, as they restore the gauge independence in 
the RG-running of neutrino masses and mixing angles. 

Next, we present analytic expressions for Majorana-neutrino self-energies in the 
gauge: 
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Likewise, we included in ( |A.11| ) the ^-dependent part of the tadpole contributions. 

Finally, for completeness, we give the charged-lepton transition amplitudes for I' —>■ I: 
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Apart from the ^-dependent part of the Higgs tadpoles included in (|A.14|) , ( |A.12| )- (|A.14| ) 



agree well with those presented in |]T2 
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